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Abstract 

I review the basics of neutrino cosmology, from the question of neutrino decouphng and the 
presence of sterile neutrinos to the effects of neutrinos on the cosmic microwave background and 
large scale structure. Particular emphasis is put on cosmological neutrino mass measurements, 



I 

o 
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ly-^ ' both the present bounds and the future prospects. 

1 Introduction 
in 

O ^ Neutrinos are the second most abundant particles in our Universe. This means they have a profound 
■ impact on many different aspects of cosmology, from the question of leptogenesis in the very early 
universe, over big bang nucleosynthesis, to late time structure formation. Here I review some general 
apects of neutrino cosmology with particular emphasis on issues relevant to cosmological bounds on the 
i neutrino mass. 

[ The absolute value of neutrino masses are very difficult to measure experimentally. On the other 
hand, mass differences between the light neutrino mass eigenstates, (mi, 7722,777,3), can be measured in 
^ ■ neutrino oscillation experiments. 

The combination of all currently available data suggests two important mass differences in the 
neutrino mass hierarchy. The solar mass difference of 6ml2 — 7.1 — 8.9 x 10~^ eV^ {3a) and the 
atmospheric mass difference Sml^ ~ 1.4 — 3.3 x 10~^ eV^ (3cr) [Tlllll2l|3] (see [HI for a recent review). 

In the simplest case where neutrino masses are hierarchical these results suggest that 7771 ~ 0, 
7772 ~ Suisoiar, and 7773 ~ ^ ^atmospheric • If the hierarchy is inverted one instead finds 7773 ~ 0, 7^72 ~ 
^'^atmospheric, and 77ii ~ ^ "^atmospheric • Howevcr, it is also possiblc that neutrino masses are degenerate, 
7r7i ~ 7712 ~ 7^73 ^ ^"^atmospheric, which casc oscillatiou experiments are not useful for determining the 
absolute mass scale. 

Experiments which rely on kinematical effects of the neutrino mass offer the strongest probe of this 
overall mass scale. Tritium decay measurements have been able to put an upper limit on the electron 
neutrino mass of 2.3 eV (95% conf.) |6j. However, cosmology at present yields an much stronger limit 
which is also based on the kinematics of neutrino mass. 

Very interestingly there is also a claim of direct detection of neutrinoless double beta decay in the 
Heidelberg-Moscow experiment [Tj |H], corresponding to an effective neutrino mass in the 0.1 — 0.9 eV 
range. If this result is confirmed then it shows that neutrino masses are almost degenerate and well 
within reach of cosmological detection in the near future. 
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Another important question which can be answered by cosmological observations is how large the 
total neutrino energy density is. Apart from the standard model prediction of three light neutrinos, 
such energy density can be either in the form of additional, sterile neutrino degrees of freedom, or a 
non-zero neutrino chemical potential. 

In section 2 I review the present cosmological data which can be used for analysis of neutrino physics. 
In section 3 I discuss neutrino physics around the epoch of neutrino decoupling at a temperature of 
roughly 1 MeV, including the relation between neutrinos and Big Bang nucleosynthesis. Section 4 
discusses neutrinos as dark matter particles, including mass constraints on light neutrinos, and sterile 
neutrino dark matter. Finally, section 5 contains a discussion. 

2 Cosmological data 

Large Scale Structure (LSS) — At present there are two large galaxy surveys of comparable size, 
the Sloan Digital Sky Survey (SDSS) CUl and the 2dFGRS (2 degree Field Galaxy Redshift Survey) 
[TTj . The SDSS is still ongoing, but will be completed very soon. It will then be significantly larger 
than the 2dF. Furthermore, a continuation in the form of a new and extended survey, SDSS-II, has 
already been approved (see h ttp: / / www.sdss.org| ) . 

Cosmic Microwave Background (CMB) — The CMB temperature fluctuations are conveniently 
described in terms of the spherical harmonics power spectrum Cf^ = (|a/mp), where ^{9,(1)) = 
Him'^imYimiO.cl)). Siucc Thomsou scattering polarizes light, there are also power spectra coming from 
the polarization. The polarization can be divided into a curl-free {{E)) and a curl {{B)) component, 
much in the same way as E and B in electrodynamics can be derived from the gradient of a scalar 
field and the curl of a vector field respectively (see for instance J2] for a very detailed treatment). The 
polarization introduced a sequence of new power spectra, but because of different parity some of them 
are exphcitly zero. Altogether there are four independent power spectra: Cf^, Cf^, Cf^, and the T-E 
cross-correlation Cf^. 

The WMAP experiment has reported data only on Cf^ and Cf^ as described in Refs. |13[ll4j. Other 
experiments, while less precise in the measurement of the temperature anisotropy and not providing full- 
sky coverage, are much more sensitive to small scale anisotropies and to CMB polarization. Particularly 
the ground based CBI [IH], DASI [HI, and ACBAR ^ experiments, as well as the BOOMERANG 
balloon experiment [IHl UHl 120] have provided useful data. 

Type la supernovae Observations of distant supernovae have been carried out on a large scale 
for about a decade. In 1998 two different projects almost simultaneously published measurements of 
about 50 distant type la supernovae, out to a redshift or about 0.8 [211 123 • These measurements were 
instrumental for the measurement of the late time expansion rate of the universe. 

Since then a, new supernovae have continuously been added to the sample, with the Riess et al. 
|2S1 "gold" data set of 157 distant supernovae being the most recent. This includes several supernovae 
measured by the Hubble Space Telescope out to a redshift of 1.7. 

Very recently, the first data has been released from the Supernova Legacy Survey (SNLS) j21], 
providing the currently largest data set of Type la supernovae. 

Other data — Apart from CMB, LSS and SNI-a data there are a number of other cosmological 
measurements of importance to neutrino cosmology. Perhaps the most important is the measurement 
of the Hubble constant by the HST Hubble Key Project, i/o = 72 ± 8 km s"^ Mpc"^ |23]. 
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3 Neutrino Decoupling 



3.1 Standard model 

In the standard model neutrinos interact via weak interactions with e"*" and e^. In the absence of 
oscillations neutrino decoupling can be followed via the Boltzmann equation for the single particle 
distribution function 1221 



df df 

where Ccou represents all elastic and inelastic interactions. In the standard model all these interactions 
are 2^2 interactions in which case the collision integral for process i can be written 



Ccoll,i(/l 



d^P2 d^Ps d^P4 



2Ei J 2£'2(27r)3 2^3(277)3 2^4(277)2 

X {2n)^6\p, +P2-P3+ P4)A(/i, /2, /s, /4)^|M|2 (2) 



where 5'|M|^2^34i is the spin-summed and averaged matrix element including the symmetry factor S = 
1/2 if there are identical particles in initial or final states. The phase-space factor is A(/i, /2, /s, f^j = 

hh{l - /l)(l - /2) - /l/2(l - /3)(1 - U). 

The matrix elements for all relevant processes can for instance be found in Ref. 28 . If Maxwell- 
Boltzmann statistics is used for all particles, and neutrinos are assumed to be in complete scattering 
equilbrium so that they can be represented by a single temperature, then the collision integral can be 
integrated to yield the average annihilation rate for a neutrino 



^U + 9l)T\ (3) 



^3 



where 



sin^ 9w + {\ + sin^ OwY for Vf. 



2,2 _ I I \2 ' ^ > (A\ 

3L + 3R-\sin^ew + {-\ + sin^ew? for z/^,. ' 

This rate can then be compared with the Hubble expansion rate 

H = IMgl/'^ (5) 

to find the decoupling temperature from the criterion H = r| j,^^^. From this one finds that T^jiy^.) ~ 
2.4 MeV, Td^u^-j-) ~ 3.7 MeV, when g^: = 10.75, as is the case in the standard model. 

This means that neutrinos decouple at a temperature which is significantly higher than the electron 
mass. When e"'"e~ annihilation occurs around T ~ me/3, the neutrino temperature is unaffected 
whereas the photon temperature is heated by a factor (11/4)^/^. The relation T^/T^ = (4/11)^/^ ~ 0.71 
holds to a precision of roughly one percent. The main correction comes from a slight heating of 
neutrinos by e+e~ annihilation, as well as finite temperature QED effects on the photon propagator 



3.2 Big Bang nucleosynthesis and the number of neutrino species 

Shortly after neutrino decoupling the weak interactions which keep neutrons and protons in statistical 
equilibrium freeze out. Again the criterion H = rL_^ can be applied to find that Tfreeze — 0.5qV^ 

' — -^freeze 

MeV 123. 
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Eventually, at a temperature of roughly 0.2 MeV deuterium starts to form, and very quickly all free 
neutrons are processed into ^He. The final helium abundance is therefore roughly given by 



+ rin/up 



(6) 

T~0.2 MeV 



rin/up is determined by its value at freeze out, roughly by the condition that n^/riplr^j-j^,^^^^ ~ 

Since the freeze-out temperature is determined by this in turn means that can be inferred 
from a measurement of the helium abundance. However, since Yp is a function of both Qi,h'^ and g^^ 
it is necessary to use other measurements to constrain Qbh"^ in order to find a bound on g^,. One 
customary method for doing this has been to use measurements of primordial deuterium to infer Qbh"^ 
and from that calculate a bound on g^,. Usually such bounds are expressed in terms of the equivalent 
number of neutrino species, N^, = p/puo, instead of g^. The exact value of the bound is quite uncertain 
because there are different and inconsistent measurements of the primordial helium abundance (see for 
instance Ref. for a discussion of this issue). The most recent analyses are jl2] where a value of 
1.7 < A^^ < 3.0 (95% C.L.) was found, which found -1.14 < AN^ < 0.73, and 03] which found that 

= 3.14lo;65 &t 68% C.L. The difference in these results can be attributed to different assumptions 
about uncertainties in the primordial helium abundance. 

Another interesting parameter which can be constrained by the same argument is the neutrino 
chemical potential, = [Xy/T jUl HHl 113 ED] • At first sight this looks like it is completely equivalent to 
constraining N^,. However, this is not true because a chemical potential for electron neutrinos directly 
influences the n — p conversion rate. Therefore the bound on from BBN alone is relatively stringent 
(—0.1 < < 1 |17]) compared to that for muon and tau neutrinos ( ^j^^^ < 7 |17j). However, as will 
be seen in the next section, neutrino oscillations have the effect of almost equilibrating the neutrino 
chemical potentials prior to BBN, completely changing this conclusion. 



3.3 The number of neutrino species - joint CMB and BBN analysis 

The BBN bound on the number of neutrino species presented in the previous section can be comple- 
mented by a similar bound from observations of the CMB and large scale structure. The CMB depends 
on N^, mainly because of the early Integrated Sachs Wolfe effect which increases fluctuation power at 
scales slightly larger than the first acoustic peak. The large scale structure spectrum depends on 
because the scale of matter-radiation equality is changed by varying N^. 

Many recent papers have used CMB, LSS, and SNI-a data to calculate bounds bounds on 
(321 Ell ESI Uni; and some of the bounds are listed in Table 13.31 Recent analyses combining BBN, 
CMB, and large scale structure data can be found in [3111121, and these results are also listed in Table 

Common for all the bounds is that = {] is ruled out by both BBN and CMB/LSS. This has the 
important consequence that the cosmological neutrino background has been positively detected, not 
only during the BBN epoch, but also much later, during structure formation. 

The most recent bound which uses the SNI-a "gold" data set, as well as the new Boomerang CMB 
data finds a limit of Ny = 4.2^\'l at 95% C.L. The bound from late-time observations is now as good as 
that from BBN, and the two derived value are mutually consistent given the systematic uncertainties 
in the primordial helium value. 

In Fig. 1 we show the currently allowed region for Qbh"^ and N^, (taken from [3^). 
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Figure 1: The 68% (dark) and 95% (light) hkehhood contours for Qhh"^ and for all available data. 
The other contours are 68% and 95% regions for BBN, assuming the ^He and D values given in 



Table 1: Various recent limits on the effective number of neutrino species, as well as the data used. 



Ref. 



Bound on 



Data used 



Crotty et al. 
Hannestad [Slj 
Pierpaoli jSSl 
Barger et al. 
Hannestad |^ 



1.4 <N^< 6.8 
0.9 <N^< 7.0 
1.9 < A^^ < 6.62 
0.9 < A^^ < 8.3 
3.0 <K< 5.9 



CMB, LSS 
CMB, LSS 
CMB, LSS 
CMB 

CMB, LSS, SNI-a 
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3.4 The effect of oscillations 



In the previous section the one-particle distribution function, /, was used to describe neutrino evolu- 
tion. However, for neutrinos the mass eigenstates are not equivalent to the flavour eigenstates because 
neutrinos are mixed. Therefore the evolution of the neutrino ensemble is not in general described by 
the three scalar functions, /j, but rather by the evolution of the neutrino density matrix, p = the 
diagonal elements of which correspond to /j. 

For three-neutrino oscillations the formalism is quite complicated. However, the difference in Ami2 
and Am23, as well as the fact that sin 2^13 ^ 1 means that the problem effectively reduces to a 2 x 2 
oscillation problem in the standard model. A detailed account of the physics of neutrino oscillations in 
the early universe is outside the scope of the present paper, however an excellent and very thorough 
review can be found in Ref. [53] 

Without oscillations it is possible to compensate a very large chemical potential for muon and/or 
tau neutrinos with a small, negative electron neutrino chemical potential jlTj. However, since neutrinos 
are almost maximally mixed a chemical potential in one flavour can be shared with other flavours, and 
the end result is that during BBN all three flavours have almost equal chemical potential. This in turn 
means that the bound on z/g applies to all species so that jHSl EZI ISHl ISH ISO] • 

le.l = ^ < 0.15 (7) 

for i = e, fi, T. 

In models where sterile neutrinos are present even more remarkable oscillation phenomena can occur. 
However, I do not discuss this possibility further, except for the possibility of sterile neutrino warm 
dark matter, and instead refer to the review |55j . 

3.5 Low reheating temperature and neutrinos 

In most models of inflation the universe enters the normal, radiation dominated epoch at a reheating 
temperature, Trh, which is of order the electroweak scale or higher. However, in principle it is possible 
that this reheating temperature is much lower, of order MeV. This possibility has been studied many 
times in the literature, and a very general bound of Trh > 1 MeV has been found |6T| l62| IHSl IM] 

This very conservative bound comes from the fact that the light element abundances produced by 
big bang nucleosynthesis disagree with observations if the universe if matter dominated during BBN. 
However, a somewhat more stringent bound can be obtained by looking at neutrino thermalization 
during reheating. If a scalar particle is responsible for reheating then direct decay to neutrinos is 
suppressed because of the necessary helicity flip. This means that if the reheating temperature is 
too low neutrinos never thermalize. If this is the case then BBN predicts the wrong light element 
abundances. However, even if the heavy particle has a significant branching ratio into neutrinos there 
are problems with BBN. The reason is that neutrinos produced in decays are born with energies which 
are much higher than thermal. If the reheating temperature is too low then a population of high energy 
neutrinos will remain and also lead to conflict with observed light element abundances. A recent analysis 
showed that in general the reheating temperature cannot be below roughly 4 MeV |65j . 

4 Neutrino Dark Matter 

Neutrinos are a source of dark matter in the present day universe simply because they contribute to 
flrn- The present temperature of massless standard model neutrinos is Tj^ o = 1.95 -ft" = 1.7 x 10~^ eV, 
and any neutrino with m ^ behaves like a standard non-relativistic dark matter particle. 
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The present contribution to the matter density of Ny neutrino species with standard weak interac- 
tions is given by 



92.5 eV 

Just from demanding that Vt^ <1 one finds the bound IHSl IHZI 



46 eV 

< — - (9) 



4.1 The Tremaine-Gunn bound 

If neutrinos are the main source of dark matter, then they must also make up most of the galactic dark 
matter. However, neutrinos can only cluster in galaxies via energy loss due to gravitational relaxation 
since they do not suffer inelastic collisions. In distribution function language this corresponds to phase 
mixing of the distribution function [HHI- By using the theorem that the phase-mixed or coarse grained 
distribution function must explicitly take values smaller than the maximum of the original distribution 
function one arrives at the condition ^ 

/CG < /i/.max = 2 (10) 

Because of this upper bound it is impossible to squeeze neutrino dark matter beyond a certain limit |68j . 
For the Milky Way this means that the neutrino mass must be larger than roughly 25 eV z/ neutrinos 
make up the dark matter. For irregular dwarf galaxies this limit increases to 100-300 eV [0111201; and 
means that standard model neutrinos cannot make up a dominant fraction of the dark matter. This 
bound is generally known as the Tremaine-Gunn bound. 

Note that this phase space argument is a purely classical argument, it is not related to the Pauli 
blocking principle for fermions (although, by using the Pauli principle /j^ < 1 one would arrive at a 
similar, but slightly weaker limit for neutrinos). In fact the Tremaine-Gunn bound works even for 
bosons if applied in a statistical sense [Oil, because even though there is no upper bound on the fine 
grained distribution function, only a very small number of particles reside at low momenta (unless there 
is a condensate). Therefore, although the exact value of the limit is model dependent, limit applies to 
any species that was once in thermal equilibrium. A notable counterexample is non-thermal axion dark 
matter which is produced directly into a condensate. 



4.2 Neutrino hot dark matter 

A much stronger upper bound on the neutrino mass than the one in Eq. Q can be derived by noticing 
that the thermal history of neutrinos is very different from that of a WIMP because the neutrino only 
becomes non-relativistic very late. 

In an inhomogeneous universe the Boltzmann equation for a coUisionless species is [TI] 

Lffl-^-^ + ^^ + ^^-0 (11) 

where r is conformal time, dr = dt/a, and = ap* is comoving momentum. The second term on the 
right-hand side has to do with the velocity of the distribution in a given spatial point and the third 
term is the cosmological momentum redshift. 

Following Ma and Bertschinger [7T] this can be rewritten as an equation for \E', the perturbed part 
of/ 

f{x\q\T) = foiq)[l + ^ix\q\T)] (12) 
In synchronous gauge that equation is 
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/c 



h + Qfi 



-fx 



fo 



C[fl 



(13) 



where = qn^ , /i = n^kj, and e = (g^ + a^m^)^/^. /c-' is the comoving wavevector. h and r] are the 
metric perturbations, defined from the perturbed space-time metric in synchronous gauge jTT] 



ds'^ = a'^{T)[—dT'^ + {6ij + hij)dx'^dx^], 



hi 



d^ke^^'"^ ( kikjh{k, r) + {kikj - -Sij)6r]{k, r) ) . 



Expanding this in Legendre polynomials one arrives at a set of hierarchy equations 



(14) 
(15) 



5 

e 

2a 
Fi 



4 2 • 

— e--h 

3 3 

^%-) 

8 . 3 , ^ 4 ■ 8 

— 9 kF^ H h + -r? 

15 15 15 5 ' 

(/F,_i - (/ + l)F,+i) 



2/ 



(16) 



For subhorizon scales {h = fj = 0) this reduces to the form 

5 = - 

e = e{--a 



2a = —9 kFs 

15 15 

Fi = ^ (/F,_i - (/ + l)F,+i) 



(17) 



One should notice the similarity between this set of equations and the evolution hierarchy for spher- 
ical Bessel functions. Indeed the exact solution to the hierarchy is 

Fiikr) jiikr) (18) 

This shows that the solution for 6 is an exponentially damped oscillation. On small scales, k > t, 
perturbations are erased. 

This in intuitively understandable in terms of free-streaming. Indeed the Bessel function solution 
comes from the fact that neutrinos are considered massless. In the limit of CDM the evolution hierarchy 
is truncated by the fact that ^ = 0, so that the CDM perturbation equation is simply 6 = —h/2. For 
massless particles the free-streaming length is A = cr which is refiected in the solution to the Boltzmann 
hierarchy. Of course the solution only applies when neutrinos are strictly massless. Once T ~ m there 
is a smooth transition to the CDM solution. Therefore the final solution can be separated into two 
parts: 1) k > t{T = m): Neutrino perturbations are exponentially damped 2) k < t{T = m): Neutrino 
perturbations follow the CDM perturbations. Calculating the free streaming wavenumber in a fiat CDM 
cosmology leads to the simple numerical relation (applicable only for T^q ^ m ^ Tq) j27j 



A 



FS 



20 Mpc /T, 



1 + log 3.9 




(19) 
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Figure 2: The transfer function T{k,t = Iq) for various different neutrino masses. The sohd (black) 
hne is for m^, = 0, the long-dashed for rriu = 0.3 eV, and the dashed for m^, = 1 eV. 

In Fig. 121 1 have plotted transfer functions for various different neutrino masses in a flat ACDM 
universe {Qm + + = !)• The parameters used were flh = 0.04, f2cDM = 0.26 — fl^, Qa = 0.7, 
h = 0.7, and n = 1. 

When measuring fluctuations it is customary to use the power spectrum, P{k,T), defined as 

P{k,T) = \6\'{T). (20) 
The power spectrum can be decomposed into a primordial part, Po{k), and a transfer function T{k, r), 

P{k,T)=Po{k)T{k,T). (21) 

The transfer function at a particular time is found by solving the Boltzmann equation for S{t). 

At scales much smaller than the free-streaming scale the present matter power spectrum is suppressed 
roughly by the factor j72| 

A^(fe) ^ AT{k,T = To) ^ 

P{k) T{k,r = To) QJ ^ ^ 

as long as Qi, ^ Qm- The numerical factor 8 is derived from a numerical solution of the Boltzmann 
equation, but the general structure of the equation is simple to understand. At scales smaller than the 
free-streaming scale the neutrino perturbations are washed out completely, leaving only perturbations 
in the non-relativistic matter (CDM and baryons). Therefore the relative suppression of power is 
proportional to the ratio of neutrino energy density to the overall matter density. Clearly the above 
relation only applies when <^ Q^n, when Q^, becomes dominant the spectrum suppression becomes 
exponential as in the pure hot dark matter model. This effect is shown for different neutrino masses in 

Fig.m 

The effect of massive neutrinos on structure formation only applies to the scales below the free- 
streaming length. For neutrinos with masses of several eV the free-streaming scale is smaller than the 
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scales which can be probed using present CMB data and therefore the power spectrum suppression 
can be seen only in large scale structure data. On the other hand, neutrinos of sub-eV mass behave 
almost like a relativistic neutrino species for CMB considerations. The main effect of a small neutrino 
mass on the CMB is that it leads to an enhanced early ISW effect. The reason is that the ratio of 
radiation to matter at recombination becomes larger because a sub-eV neutrino is still relativistic or 
semi-relativistic at recombination. With the WMAP data alone it is very difficult to constrain the 
neutrino mass, and to achieve a constraint which is competitive with current experimental bounds it is 
necessary to include LSS data from 2dF or SDSS. When this is done the bound becomes very strong, 
somewhere in the range of 1-1.5 eV for the sum of neutrino masses, depending on assumptions about 

priors [HEIllHailDllHIllZaiZilZZllZHllZilHDl 

The bound can be strengthened even further by including data from the Lyman-a forest, measured 
at large redshifts. This was done for instance in [7^ I77j . where bounds as strong as 0.4 eV were derived. 
However, the systematic errors in extracting the matter power spectrum from the Lyman-a flux power 
spectrum are at present not understood at a level where this bound can be claimed to be robust. 

4.2.1 Parameter degeneracy with the dark energy equation of state 

One caveat of most cosmological parameter analyses used to probe particle physics is that they are done 
with the minimal cosmological standard model. In principle there could easily be additional parameters 
which are important. One such example was described in [HU], where it was shown that there is an 
almost perfect degeneracy between the neutrino mass and the equation of state of the dark energy, w. 
This degeneracy is shown in Fig. El An increasing ''^u can be compensated by decreasing w. While 
for low neutrino masses a cosmological constant {w = —1) is allowed, for high neutrino masses only 
dark energy models in the phantom regime {w < —1) are allowed. 



1 , , , \ , , , , f 




Sm^, (eV) 

Figure 3: 68% and 95% allowed contours as a function of neutrino mass and dark energy equation of 
state using WMAP, SDSS, HST, and SNI-a data. 

The reason for the degeneracy is that when Q^, is increased, Qm must be increased correspondingly in 
order to produce the same power spectrum. However, when w = —1 sua increasing Qm quickly becomes 
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incompatible with the supernova data. This can be remedied by simultaneously decreasing w because 
of the well-known Q^, w degeneracy in the supernova data. For the particular data set used in this case 
the 95% C.L. bound changes from 0.65 eV with w = —1 to 1.5 eV with free w. 

4.2.2 Combining measurements of niy and A^^. 

The limits on neutrino masses discussed above apply only for neutrinos within the standard model, i.e. 
three light neutrinos with degenerate masses (if the sum is close to the upper bound). However, if there 
are additional neutrino species sharing the mass, or neutrinos have significant chemical potentials this 
bound is changed. Models with massive neutrinos have suppressed power at small scale, with suppression 
proportional to VLy/VLm- Adding relativistic energy further suppresses power at scales smaller than the 
horizon at matter-radiation equality. For the same matter density such a model would therefore be 
even more incompatible with data. However, if the matter density is increased together with and 
Ny, excellent fits to the data can be obtained. This effect is shown in Fig. 0] 



1000 
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Figure 4: Power spectra for ACDM models with Vti, = 0.05, f2 
and a common large-scale normalization. The full line is for 

= 0.05, fi„ = 0.25, Ny = 3, dotted is for = 0.05, fi„ 

= 0.05, = 0.35, A^^ = 8 (from \Ul). 

The effect on likelihood contours for {Q,y, N^) can be seen in Fig. which is for the case where A^i, 
species the total mass equally. 

A thorough discussion of these models can be found in Refs. |7H ITS]. 

4.2.3 Future neutrino mass measurements 

The present bound on the sum of neutrino masses is still much larger than the mass difference, Am23 ~ 
0.05 eV [HHlIll; measured by atmospheric neutrino observatories and K2K . This means that if the sum 
of neutrino masses is anywhere close to saturating the bound then neutrino masses must the almost 



--l,h = 0.7, ris = 1, and A^,.,massive = 1 
= 0, flm = 0.25, Nu = 3, dashed is for 
0.25, Ny = 8, and long-dashed is for 
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Figure 5: Likelihood contours (68% and 95%) for the case of N^, neutrinos with equal masses, calculated 
from WMAP and 2dF data (from [71]). 

degenerate. The question is whether in the future it will be possible to measure masses which are of 
the order Am23, i.e. whether it can determined if neutrino masses are hierarchical. 

By combining future CMB data from the Planck satellite with a galaxy survey like the SDSS it 
has been estimated that neutrino masses as low as about 0.1 eV can be detected j83| IMj . Another 
possibility is to use weak lensing of the CMB as a probe of neutrino mass. In this case it seems likely 
that a sensitivity below 0.1 eV can also be reached with CMB alone jHO]- Possibly the sensitivity could 
be increased even further with future galaxy cluster surveys PT| . 

As noted in Ref. |H1] the exact value of the sensitivity at this level depends both on whether the 
hierarchy is normal or inverted, and the exact value of the mass splittings. 

4.3 Neutrino warm dark matter 

While CDM is defined as consisting of non-interacting particles which have essentially no free-streaming 
on any astronomically relevant scale, and HDM is defined by consisting of particles which become non- 
relativistic around matter radiation equality or later, warm dark matter is an intermediate. One of the 
simplest production mechanisms for warm dark matter is active-sterile neutrino oscillations in the early 
universe CDll 1103 1103 IMl • 

One possible benefit of warm dark matter is that it does have some free-streaming so that structure 
formation is suppressed on very small scales. This has been proposed as an explanation for the apparent 
discrepancy between observations of galaxies and numerical CDM structure formation simulations. In 
general simulations produce galaxy halos which have very steep inner density profiles p oc r", where a ~ 
1 — 1.5, and numerous subhalos [021 EH!- Neither of these characteristics are seen in observations and the 
explanation for this discrepancy remains an open question. If dark matter is warm instead of cold, with 
a free-streaming scale comparable to the size of a typical galaxy subhalo then the amount of substructure 
is suppressed, and possibly the central density profile is also flattened [511 jSHl HZl ISHl ISSl HOOj . In 
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both cases the mass of the dark matter particle should be around 1 keV jlU8t llUyj . assuming that it is 
thermally produced in the early universe. 

On the other hand, from measurements of the Lyman-a forest flux power spectrum it has been 
possible to reconstruct the matter power spectrum on relatively small scales at high redshift. This 
spectrum does not show any evidence for suppression at sub-galaxy scales and has been used to put a 
lower bound on the mass of warm dark matter particles of roughly 500 eV An even more severe 

problem lies in the fact that star formation occurs relatively late in warm dark matter models because 
small scale structure is suppressed. This may be in conflict with the low-Z CMB temperature-polarization 
cross correlation measurement by WMAP which indicates very early reionization and therefore also early 
star formation. One recent investigation of this found warm dark matter to be inconsistent with WMAP 
for masses as high as 10 keV jMj . 

The case for warm dark matter therefore seems quite marginal, although at present it is not defini- 
tively ruled out by any observations. 



5 Discussion 

Here I have reviewed some of the basics of neutrino cosmology with particular emphasis on the role 
that neutrinos play in cosmological structure formation. Exactly because neutrinos play such a crucial 
role in cosmology it is possible to use cosmological observations to probe fundamental properties of 
neutrinos which are not easily accessible in lab experiments. Particularly the measurement of absolute 
neutrino masses from CMB and large scale structure data has received significant attention over the 
past few years. 

Another cornerstone of neutrino cosmology is the measurement of the total energy density in non- 
electromagnetically interacting particles. For many years Big Bang nucleosynthesis was the only probe 
of relativistic energy density, but with the advent of precision CMB and LSS data it has been possible 
to complement the BBN measurement. At present the cosmic neutrino background is seen in both 
BBN, CMB and LSS data at high significance. 

Finally, cosmology can also be used to probe the possibility of neutrino warm dark matter, which 
could be produced by active-sterile neutrino oscillations. 

In the coming years the steady stream of new observational data will continue, and the cosmological 
bounds on neutrino will improve accordingly. For instance, it has been estimated that with data from 
the upcoming Planck satellite it could be possible to measure neutrino masses as low as 0.1 eV, and 
possibly even lower than that if data from future cluster surveys is used. 

Neutrino cosmology is certainly an exiting field which is destined to play a key role in neutrino 
physics in the coming years. 
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